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Abstract 

We have computed the loop-induced processes of neutrahno annihilation into two photons and, for 
the first time, into a photon and a Z° boson in the framework of the NMSSM. The photons produced 
from these radiative modes are monochromatic and possess a clear "smoking gun" experimental signa- 
ture. This numerical analysis has been done with the help of the SloopS code, initially developed for 
automatic one-loop calculation in the MSSM. We have computed the rates for different benchmark 
points coming from SUGRA and GMSB soft SUSY breaking scenarios and compared them with the 
MSSM. We comment on how this signal can be enhanced, with respect to the MSSM, especially in 
the low mass region of the neutralino. We also discuss the possibility of this observable to constrain 
the NMSSM parameter space, taking into account the latest limits from the FERMI collaboration on 
these two modes. 



Introduction 



The existence of cold dark matter (CDM) is supported by many geometrical and dynamical observations 
coming from cosmology and astrophysics, but its detection needs to be confirmed. Its properties, such as 
mass, spin, and interactions are still to be determined. A long term experimental effort has been devoted 
to its detection, from direct methods, where a dark matter particle is expected to impinge on a nuclei 
in a target material, or more generally from the products of its self-annihilation in outer space (indirect 
detection) . 

Recently the CDMS [Ij, CoGeNT [2l|3] experiments and very recently CRESST |4j have pointed out results 
whose explanation could be interpreted as hints in favour of the existence of low mass dark matter par- 
ticles, although the recent result of the XENONIOO disfavors this possibility [5] and the CDMS- 1 1 results 
seems to be in contention with CoGeNT [IJ. The interesting piece of information is that these results could 
be reconciled with the long standing claims of the DAMA/LIBRA experiment that dark matter was detected 
through its annual modulation [6j. However these potential signals of dark matter are still an open 
debate and no firm claim of dark matter signal can be established, since there is no general agreement 
between the different experiments. If a dark matter hypothesis is to be put forward to interpret these 
signals, the studies points to a low mass WIMP. Thus some attempts were made to account for them, 
for example with neutralinos lighter than, say, 15 GeV. In recent publications the authors of [7-,12j con- 
sidered and investigated the occurrence of light neutralino within the Minimal Supersymmetric Standard 
Model (MSSM) and Next-to- Minimal Supersymmetric Standard Model (NMSSM) while respecting several 
constraints. It was shown that an MSSM explanation for the signal is seriously challenged by collider 
searches and measurements in the flavour sector |i71l8l[T3l[T^ . 

The Next-to-Minimal Supersymmetric Standard Model is a well-motivated extension of the MSSM con- 
structed by adding a singlet superfield S in the Higgs sector. This has been advocated to provide an 
elegant solution to the so-called /i problem of the MSSM but it also renders the Higgs "little fine tuning 
problem" less severe |15) . Moreover it leads to a richer spectrum than the MSSM : two additional Higgs 
bosons (one CP-even and one CP-odd neutral Higgs) remain after ElectroWeak Symmetry Breaking 
(EWSB) and the neutralino sector is enlarged by a fifth neutralino. This last possibility has important 
consequences concerning dark matter studies, that can differ significantly from the MSSM, especially 
when the lightest neutralino composition is dominated by the singlet component, the singlino. Moreover 
the occurrence of very-light singlet pseudoscalar in the Higgs spectrum and at the same time evading 
the LEP bounds is still possible thanks to reduced couplings to standard model particles. These facts 
open new possibilities in relic density calculation as well as direct detection prospects as opposed to the 
MSSM. In particular the thermally produced neutralinos can satisfy the measured CDM present relic 
abundance while being very light. This implies that the pair annihilation rate of thermal relics roughly 
scales as \/m^. All these observations make a light supersymmetric candidate for dark matter more 
plausible in the NMSSM than the MSSM. It is then appealing to investigate the self annihilation of 
dark matter in our galaxy since one expects light neutralinos to give significantly enhanced rates for the 
indirect detection signals compared to the standard case. 

Numerous studies have been devoted to relic density and direct detection phenomenology in the NMSSM, 
few work has been dedicated to indirect detection of dark matter (annihilation of a pair of dark matter 
particle) through the "direct" annihilation into primary monochromatic photons ^16j . These primaries 
photons have the advantage of being less affected by astrophysical uncertainties over other kind of mes- 
sengers. Moreover their spectrum would reveal a sharp peak at an energy ~ My- corresponding to the 
mass of the dark matter particle for the 77 final state and ~ M^(l — M^/4M^) for ^Z^ , since in the 
galactic halo f/c ~ 10~^. Provided one has a good enough detector energy resolution and sensitivity, the 
flux from these primaries photons will be clearly distinctive from the astrophysical background or diffuse 
emission. However modeling the dark matter halo is still needed since the number density of dark matter 
particles enters the calculation. The FERMI collaboration did not report any line observation and has 
released upper limits on the direct annihilation (av)^^ and {ctv)^z [TZj instead. The recently installed 
AMS [18| detector on the International Space Station may also shed light on these channels. 
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The complete computation of the loop-induced annihilation into photons has already been performed 
in |19fl22) for the MSSM. In |16) the one-loop amplitudes for NMSSM neutralino pair annihilation into 
two photons and two gluons have been given (adapted to the NMSSM case from the formulas given in |19] ) 
and prospects for the indirect detection of the monochromatic gamma-ray line were also investigated. 
In light of the latest limits of the FERMI collaboration on monochromatic gamma ray signal and the 
recent activity in the dark matter community, we propose to revisit the two gammas mode and provide 
for the first time results on the one. We investigated to which extent the NMSSM rates could 
differ from the MSSM and for the numerical study we focused on mSUGRA and GMSB benchmarks 
points given in |23p24). This is a case study if one wants to discriminate the NMSSM from the MSSM 
with dark matter related observables, which basically means to which amount the neutralino sector is 
sensitive to the extended Higgs sector. We further explored in which case, specific to the NMSSM, the 
spectral line can be enhanced and if such mechanisms can be constrained by astrophysical measurements. 
Constraining these mechanisms is very interesting since it also impacts the relic density calculation and 
direct detection predictions. 

These rates were computed with the help of the SloopS program, an automatic code for one-loop calcu- 
lation in the context of the SM and the MSSM |25|I26|. This code has already been used for accurate relic 
density predictions at next-to-leading order in the MSSM |27fl29| and also to the numerical computation 
of the indirect channels XiX? ~^ 11 il'^^ ^iid two gluons mode |22) . Therefore this work is also a good 
exercise to test the implementation of the NMSSM in SloopS. On the technical level, as the dark matter 
particles are moving at relatively small velocity, the Gram determinant, which is a key ingredient for 
calculating the loop integrals, vanishes and this results in numerical instabilities. This has been handled 
with a procedure called segmentation [22]. We will come back to the numerical and technical details 
of the implementation later. The outline of the paper is as follows, in section one we quickly review 
the NMSSM model and its parameter space, in section two we describe the SloopS code and the im- 
plementation of the NMSSM. A particular attention will be paid to the implementation of non-linear 
gauge-fixing in the NMSSM. The third section will be devoted to scrutinise the additional contributions 
brought when going from the MSSM to the NMSSM. In the fourth section we will give the rate for the 
77 and ^Z^ channels for some SUGRA and GMSB scenarios and compare it to an equivalent MSSM 
spectrum. In the following section we will discuss on some possible ways of increasing the signal, specific 
to the NMSSM, and discuss the relevance of using the sharp gamma lines as a constraining observable, 
taking into account the latest published limits. Finally we will draw our conclusions. 

1 Overview of the NMSSM 

In the NMSSM the Higgs term of the superpotential involving the two Higgs doublet is modified and a 
singlet term is addecj^] [15j, 

Wf^MSSM = Wl=^sM + ^SHuiid + (1.1) 

The MSSM \x bilinear term is now absent from the superpotential and has been replaced by the trilinear 
coupling of the singlet with the two Higgs doublets. The VEV of the singlet generates an effective \i 
parameter with respect to the MSSM, which is then naturally of order the EW scale [15| . 

/ieff = As (1.2) 

where s = {S) is the VEV of the Higgs singlet. The soft-SUSY breaking Lagrangian is also modified 
according to 

+ (XAxHu-HdS+^KA^S^ + h.c) (1.3) 

^We stick to the "Zs-invariant NMSSM", where any dimensionful parameters in the superpotential are forbidden. 



2 



Given Mz and using conditions coming from the minimisation of the Higgs potential, one can choose six 
independent parameters for the Higgs sector 



A, K,Ax,A^,fJ,eS,ti3 



[lA) 



where = tan/3 = Vu/vd, the ratio of the two Higgs doublet VEV's : (H^) = v^, {H^) = Vii . This 
is in contrast with the MSSM where two parameters are needed, e.g, tan/3 and the mass of the Higgs 
pseudoscalar M^o, once the requirement of vanishing tadpoles has been imposed. After EWSB the 
NMSSM Higgs sector contains three neutral scalar fields. Hi, H2, H-^ and two pseudoscalar neutral ones, 
Ai , A2 as well as a charged Higgs . In the neutralino sector the additional singlino mixes with the 
bino, wino and Higgsinos fields. The neutralino mass matrix is therefore a 5 x 5 one which is diagonalised 
with a unitary matrix N. In the basis = {—iB, —iW^, H^, H^, S^) the neutralino mass matrix reads. 



V 



Ml 


-CfsswMz 
spswMz 






M2 
c/^cwMz 
-spcwMz 




-CfsswMz 
cpcwMz 


— /^eff 
-\Vu 



spswMz 

-Si3CwMz 





\ 


-\Vu 

2ks / 



(1.5) 



The lightest supersymmetric particle (LSP) can then be expressed as a linear combination of the five 
gauge eigenstates 

X? = A^ii^ + N12W + Ni^Hi + NiiH2 + N15S (1.6) 



For a pure state the singlino mass is 



rrii 



2ks 



(1.7) 



The parameter space of the NMSSM can then be described, in addition to the six Higgs sector parameters, 
by the same as the MSSM, namely the soft masses for sfermions Mj, trilinear couplings Af and gaugino 
masses Mi^2,3 which are respectively the U{1), SU{2) and SU{3) soft parameters. If the LSP has a 
dominant singlino component, it can efficiently annihilate through light singlet Higgses as well as light 
pseudoscalar Higgs singlet |16 | l30 j . The MSSM limit is recovered when A, k — >• and the equivalent of 
the mass (squared) of the only physical CP-odd scalar A^ of the MSSM is given by 



M 



AO 



2Xs{Ax + Ks) 
sin2/3 



(1.8) 



This is an important proviso when we will compare the NMSSM result with an "equivalent" MSSM 
spectrum. 



2 Set-up of the automatic calculation 

One loop processes calculated via the diagrammatic Feynman approach imply the calculation of hundreds 
of Feynman diagrams and a hand calculation is proned to numerous errors. A high-level of automation 
is therefore highly desirable, especially if one wants to build a general purpose code. The SloopS code 
|22 | I25 | [25] has been developed in this purpose and applied to astrophysics |22| . cosmology |27H29| and 
also collider physics |I26|I31). both in the SM and MSSM. The implementation of the NMSSM has been 
carried out with LanHEP [32] which generates the complete set of NMSSM vertices once the Lagrangian 
is specified. The calculation of the process XiXi ~^ requires the field renormalisation 6Z^^^ which is 
generated from the (tree-level) XiXi^ vertex through a Z — j one-loop transition. This renormalisation 
constant is generated in LanHEP by shifting the appropriate fields and defined in the on-shell scheme. 
The output files are then written in the FormCalc |33| conventions which handle the calculation of 
the cross section. One nice feature of SloopS is the use of a generalised non- linear gauge fixing |34j, 
adapted to the supersymmetric case |22|I25]. and for this particular work we have extended it to the 



3 



specific NMSSM Higgs sector. Compared to the MSSM, the gauge-fixing Lagrangian wih depend on two 
additional parameters since the Higgs sector is enlarged with one more scalar and one more pseudoscalar. 
This allows us to perform a non-trivial check of the gauge independence of the result through the variation 
of eleven non-linear gauge parameters. The gauge fixing Lagrangian reads in a general form 

where the non-linear functions of the fields F are given by 
F+ = {d^-ieaA^-igcwpZ^W^' + 

+iiwi^ {v + 6iHi + 62H2 + SsHs + i{kG° + piAi + ^2^2)^ G+ 
F^ = d^Z^' + iz^[v + hHi + ~e2H2 + hH^G'' 



(2.2) 
(2.3) 



= d^.A^' (2.4) 

where and G^ are respectively the neutral and charged goldstones. The parameters a, /3 • • • £3 are 
generalised gauge fixing parameters. 

The ghost Lagrangian is derived by requiring that the full effective Lagrangian is invariant under 
BRST transformations. This implies that the full quantum Lagrangian, with Cc and Cch respectively 
the classical and ghost Lagrangians, 

Cq = Cc + Cgf + Cch (2.5) 

be such that (5brs>Cq = and hence (5brs^gf = — (^brs^Cg/i |35i |. The BRST transformation for the gauge 
fields can be found for example in |35| . The NMSSM specific transformations for the scalar fields can be 
found in the Appendix. Within this particular gauge fixing we can set ^w,z,A = 1 (avoiding complicated 
tensor structure for the gauge bosons propagators) and keep the possibility to check the gauge invariance 
of the result, at the expense of adding new vertices to the model. 

The version of LoopTools |36| we used is a modified one which tackles the problem of inverse Gram 
determinant. This kinematic quantity vanishes when the velocity v is equal to zero. Consequently the 
reduction algorithm expressing the tensorial loop integrals onto a basis of scalar ones breaks down. The 
details of this procedure can be found in |22j. In short this method uses the particular kinematics at 
f = to reduce, for example, scalar box integrals to a sum of triangle ones. The neutralino mass matrix is 
diagonalised numerically with a complex unitary matrix N giving positive eigenvalues for the neutralino 
masses. 



3 Additional NMSSM contributions 

In the MSSM, depending on the nature of the neutralino, the predicted rates for (crw)^^/^ can be quite 
different. The highest ones are reached when the LSP is mostly wino or higgsino like. In the former case 
the dominant contributions are from loops with gauge bosons, since in the wino case the most important 
coupling is XiXiW^- For a Hig gsino Xj* the dominant contribution would be vertex corrections through 
a boson in the s-channel as the XiXiZ coupling is proportional to 

g^o^oz ^ Nfs - Nf^ (3.1) 

The third possible case, a bino-like neutralino, would annihilate mainly through box diagrams containing 
right-handed sfermions since the bino, as the superpartner of the B boson, couples to sparticles possessing 
the highest hypercharge. A bino like LSP gives generically lower rates than the two other cases as its 
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Figure 1: Additional NMSSM diagrams with s-channel pseudoscalar exchange for the processes XiXi ~^ 
jj{Z^). The label f stands for a SM fermion and to a chargino. 



couplings to other particles are weaker. Therefore if one wants to stick to a mostly bino neutralino 
while giving good prospects for the spectral photon lines, significant mixing with the other components 
is needed. 

In going from the MSSM calculation to the NMSSM one, two major differences appear : on the one hand 
the NMSSM couplings now depend on the parameters A and k, and on the other hand new diagrams are 
to be computed. Concerning the case of the modified couplings, if the NMSSM neutralino is a mostly 
bino, wino or higgsino, we do not expect major differences with the MSSM, as long as the singlino 
component is negligible. As the neutralino is a Majorana fermion, at vanishing relative velocity the 
LSP pair carries pseudoscalar quantum numbers. Therefore the NMSSM additional diagrams are the 
ones with pseudoscalar Ai,A2 s-channel exchange. The relevant diagrams are depicted in Fig. [TJ Note 
that the diagram with the chargino loop should be understood as containing both "flavours" of charginos 
since the couplings Ai^xfxJ ^re non-diagonal. Note also that the couplings Ai^fifj are also non-diagonal 
in family space because we stick to real soft SUSY breaking terms. Hence one could think that the 
NMSSM brings also new contributions containing loops with different flavours of sfermions. In fact their 
amplitudes cancel when summing over all the flavour^. Furthermore due to our particular gauge-fixing 
Lagrangian in Eq. ()2.ip , we generate couplings of the pseudoscalars to the charged goldstones and to 
the charged ghosts c^,c^, proportional to a and pi^2- These additional couplings are not present in the 
linear gauge. However, just like the case of the sfermions coupling to Ai^2 the related amplitudes vanish 
identically. This is however not surprising since there exists no "tree-level" coupling of the pseudoscalars 
to a pair of gauge bosons. In turn, the process with 77 will involve only the non-linear gauge parameter 
a and the process with ^Z^ the parameter /3 in addition. Obviously the final result should not depend 
on these two parameters, they are used as a check on the gauge invariance of the cross sections, which 
we ascertained. 

The use of the non-linear gauge fixing Lagrangian of Eq. ()2.ip in the 't Hooft-Feynman gauge [^a = 
iz = = 1) cilso enables us to cancel, for example, the W^G^^ vertex, by setting a = —1, which 
cancels a corresponding piece of the original trilinear sector of the Lagrangian, leading to a vanishing 
total W~^G~'y trilinear coupling (see [37] and references therein). We can perform the same trick for the 
coupling W~^G~ by setting /3 to an appropriate value. Finally, as a further check on the implementation 
we took the input parameters for the six scenarios presented in [22], compared with our results by taking 
the MSSM limit, and found a perfect agreement. 

The FERMI collaboration has already provided their limits |17) on the 77 and jZ^ channels for different 
photon energies E^, hence we did not model the propagation of the photon signal and we restricted 
ourselves to the discussion of {av)^^ and {ov).^z- 

^This is of course also the case in the MSSM. 
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4 Benchmark points and comparison with the MSSM 



As a first appfication of our code we computed the gamma-ray lines for the SUGRA and GMSB benchmark 
points from [ 23p24) and NMSSMTools_2 .3.5 |38) . For a discussion about these benchmarks we refer to their 
respective pubhcation ^23»24j . We then compared with the present hmit given by FERMI and against 
the MSSM calculation, performed with the SloopS code with a MSSM model file, like in [22], to see how 
much precision one would need to distinguish them. The input parameters for the MSSM parameters are 
the same as in the NMSSM taken in the decoupling limit A, k — )• and using Eq. p.8p . It is important 
to note that, in order to preserve gauge invariance, we only took the input parameters obtained after 
the RGE running (therefore DR parameters) to the low-energy scale and considered them as physical 
parameters, no mass correction or effective couplings have been considered. The self-annihilation of the 
neutralino is mainly driven by its mass and composition, so the first question to ask when comparing the 
MSSM with the NMSSM is whether their mass and composition are similar. 



4.1 SUGRA 

The SUGRA input parameters are taken from the version 2.3.5 of NMSSMTools as they have been 
substantially modified since the publication of |23) . The most relevant parameters are given in Table [TJ 
The first three scenarios (PI, P2, P3) exhibit a mostly bino-like neutralino and are quite similar, as far 



Parameter 


SUGRA-Pl 


SUGRA-P2 


SUGRA-P3 


SUGRA-P4 


Ml 
Ma 

MjR 

'1,2 


211.62 
391.86 
968.62 
10.000 
257.51 


211.62 
391.86 
968.62 
10.000 
257.51 


211.86 
392.51 
938.41 
10.000 
256.38 


333.80 
613.91 
-206.17 
3.0413 
702.15 




221.13 


221.13 


221.21 


698.94 




986.12 


986.12 


985.64 


1674.5 


565.23 


565.23 


590.90 


527.83 




981.60 


981.60 


981.17 


1609.0 


966.15 


966.15 


966.20 


1606.6 


A, 

Ar 

A^ 

At 


-1777.5 
-1763.0 
-2611.8 
-1431.5 


-1777.5 
-1763.0 
-2611.8 
-1431.5 


-1698.7 
-1682.8 
-2535.2 
-1362.6 


-2531.2 
-2524.8 
-3730.5 
-1765.8 



MSSM specific parameter 



M 



A" 



943.65 



943.65 



935.19 



702.11 



NMSSM specific parameters 



A 


0.1 


0.1 


0.4 


0.49 


K 


0.1089 


0.1089 


0.3037 


0.0554 


Ax 


-963.91 


-963.91 


-620.23 


-684.51 


A. 


-1.5893 


-1.0934 


-11.158 


151.73 



Table 1: SUSY input parameters for the SUGRA benchmarks. Masses are in GeV. 



as only the neutralino/chargino sector is concerned. In the four cases the lightest neutralino is heavy 
enough (?n^o ~ 211 GeV) such that the ^Z^ final state is open. As the first three benchmark scenarios 
give a bino-like LSP, we can already anticipate that the NMSSM and MSSM will give similar results 
as they have a very similar composition. However we expect that the fourth scenario P4 will give a 
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clearly distinctive signature since in this case the hghtest NMSSM neutralino x? is 99% sin fflino with 
mass m^o ~ 60 GeV. The final results can be found in Table [2j 



Model 


(GeV) 


(av)^^ X 10^'^[cm^s-^] 


{av)^z X 10^'^[cm^s-^] 


SUGRA-P1/P2 


NMSSM 


210.74 


3.9840 


1.7391 


MSSM 


210.74 


3.9836 


1.7393 


SUGRA-P3 


NMSSM 


210.92 


4.0103 


1.7326 


MSSM 


210.93 


4.0074 


1.7348 


SUGRA-P4 


NMSSM 


59.963 


2.715710"^ 


2.101210"^ 


MSSM 


200.43 


61.866 


205.28 



Table 2: Rates for the loop-induced annihilation of two neutralinos into 77/^*^ for the SUGRA bench- 
marks. 

We see that, unless the NMSSM neutralino has a radically different composition from the MSSM case, 
these two frameworks are almost indistinguishable. Needless to say that this is also a good indication 
of the implementation of the NMSSM model since in the decoupling limit (which is the case for these 
two scenarios) we should recover the MSSM. The benchmarks SUGRA-Pl and SUGRA-P2 give equal 
results because they only differ by the value of A^^, which has an impact only on the mass of the lightest 
pseudoscalar, which is not relevant in these cases. A difference is clearly visible when the NMSSM 
neutralino is mostly singlino, see the SUGRA-P4 case, whereas the MSSM neutralino is higgsino-like. 
Obviously for this latter case the comparison is not very meaningful from an experimental point of view 
since the inputs would be the energy of the photon and the rate (fiu), which are clearly different 
in this benchmark point. Considering only the first three NMSSM SUGRA scenarios, the gamma-line 
observable does not permit us to discriminate them. The distinct features of the NMSSM Higgs sector are 
quite decoupled from the rest because A and k are quite small for the first two scenarios. The situation is 
different for the third scenario where A and k take higher values and the NMSSM Higgs sector exhibits a 
different spectrum than the MSSM one because, in this case, the singlet component mixes notably with 
the doublets. Nevertheless, we recall that as the neutralino is a Majorana particle, it forms a pseudoscalar 
state such that its couplings to scalar Higgses are suppressed at vanishing relative velocities, and therefore 
the LSP is mainly sensitive to the pseudoscalar part of the Higgs spectrum. Moreover the singlino-like 
neutralino is very heavy compared to the rest of the neutralino spectra, such that it does not bring 
significant mixing. Thus the couplings of the lightest neutralino does not depart much with respect to 
the MSSM case. In addition a MSSM bino-like neutralino does not couple much to Higgs at low t^, such 
that in the three bino scenarios considered the annihilations are driven by box diagrams with sfermions, 
in particular staus. Let us now comment the SUGRA-P4 scenario. The NMSSM benchmark point gives 
a LSP which is almost a pure singlino and couples very feebly to other particles. Hence we obtain a very 
small rate, compared to the Higgsino-like LSP in the P4 MSSM benchmark, since in this case the pair 
annihilation is quite efficient as the LSP has SU{2) quantum numbers. 

All in all, taking into account the latest limits provided by FERMI [iTJ, none of these benchmarks points 
are excluded and far away from FERMI sensitivity. As these benchmarks are representative of the NMSSM 
phenomenology this means and confirms the FERMI collaboration statement that the present limits on the 
gamma-ray line are too weak compared to the typical cross sections for a conventional thermal NMSSM 
LSP. We also remark that the situation is worse for a mostly singlino neutralino. This is quite generic for 
a singlino Xi unless some specific mechanism is at play, we will come back to this point later. As a final 
comment to this section we also observe that the 7Z mode gives generically smaller rates (in both the 

^The input value of t/3 is equal to 10 at the Mz scale 
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MSSM and NMSSM) than the 77 one, as usuahy claimed, except for the SUGRA-P4 MSSM scenario. 



4.2 GMSB 

We now turn to the GMSB benchmark points. These were taken once again from NMSSMTools_2 . 3 . 5 
package |38] . The most relevant ones are depicted in Table |3l The original GMSB phenomenology of the 
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vjiViO-D-r Z 


vjiVio-D-_r 


PA/TCIR P4 
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1 "id fiO 

ioD.oy 
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00^.0 1 


1 '^fi'^ n 

lOUO.U 




8.4821 


1.6277 


1.5982 


1.9000 


50.055 


Mm 
'1,2 


692.31 


689.42 


138.22 


133.34 


622.57 




686.00 


689.20 


138.17 


133.28 


429.74 


M~R 

"1,2 


231.81 


231.70 


771.61 


777.47 


2439.5 




187.05 


176.49 


631.97 


660.97 


20361.8 




228.40 


228.27 


767.45 


773.63 


2413.4 




227.61 


228.23 


767.36 


773.52 


2238.0 


A, 


-424.24 


-321.01 


-54.511 


-46.704 


-208.15 




-423.01 


-318.15 


-54.333 


-46.573 


-184.48 


A, 


-2115.1 


-1990.8 


-439.67 


-398.57 


-166.66 


At 


-1558.5 


-1265.4 


-328.32 


-314.23 


-143.41 




VISSM specific parameter 


M40 


1756.1 


2884.8 


875.97 


714.730 


1061.3 


NMSSM specific parameters 


A 


0.002 


0.5 


0.6 


0.6 


0.01 




0.0045 


0.4351 


0.4540 


0.3969 


-0.0007 


Ax 


-56.263 


-449.46 


23.012 


13.083 


114.94 


Ak 


-158.89 


-2278.7 


0.0315 


0.7354 


0.0048 



Table 3: SUSY input parameters for the GMSB benchmarks. Masses are in GeV. 



NMSSM has been studied in |2^. The spectral lines results can be found in Table SI The GMSB-Pl and 
P2 benchmarks give both almost similar results with m^o ~ 472 GeV. We observe no major discrepancy 
between the NMSSM and MSSM results, even in scenario P2 where the NMSSM sector is not in the 
decoupling regime since A and n take quite high values. Nevertheless it still has no impact because the 
LSP is 99% bino-like and has MSSM-like couplings. Note that for these two points the LSP mass is 
outside the energy range published in |17|I39). For benchmarks P3 and P4 the differences between the 
NMSSM and the MSSM are tiny, already at the level of the mass of the lightest neutralino m^o , in turn 
the rates are slightly different. This is mostly due to the fact that the signal is computed at a different 
center of mass energy. Even though it would correspond to different signal/energy bins, we can anyway 
compare them since the mass difference is quite small and require a high level of precision. Therefore even 
discriminating between the two models is also extremely challenging, both in 77 and '^Z^ final states. 
As well as the mSUGRA benchmark points where the LSP is mostly bino-like, the annihilation is driven 
by box diagrams involving right-handed sfermions, and in particular the right-handed sleptons since they 
possess the lower masses of the sfermion spectrum. The benchmarks P3 and P4 give an enhanced rate 
because the right-handed slepton masses are close to the mass of the neutralino. In addition the bino-like 
neutralino in the P4 scenario has a substantial higgsino component, increasing its couplings to the 
boson. Finally the scenario P5 gives an extremely tiny signal where the LSP is mostly singlino and the 
s-channel exchange of the singlet pseudoscalar is extremely suppressed since its mass is around 1 GeV. 
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Model 


(GeV) 


{av)^^ X 10^'J[cnr^s-i] 


{(jv)^z X 10^'^[cm^s-i] 


GMSB-Pl 


NMSSM 


472.42 


0.2915 


0.1176 


MSSM 


472.42 


0.2915 


0.1176 


GMSB-P2 


NMSSM 


472.47 


0.2981 


0.1185 


MSSM 


472.47 


0.2981 


0.1185 


GMSB-P3 


NMSSM 


133.01 


12.862 


6.0282 


MSSM 


133.02 


12.827 


6.0466 


GMSB-P4 


NMSSM 


132.50 


14.490 


7.1262 


MSSM 


132.55 


14.460 


7.1544 


GMSB-P5 


NMSSM 


196.84 


2.626410-1" 


5.559610-11 


MSSM 


496.83 


1.0884 


0.5055 



Table 4: Rates for the loop-induced annihilation of two neutralinos into for the GMSB benchmarks. 



5 Spectral lines from light neutralino annihilation 

We have seen in the previous section that a mostly bino LSP in the NMSSM mimics the MSSM result and 
gives quite low predictions with respect to the present experimental limits. We propose in this section 
to investigate how the signal can be improved in the NMSSM, with features distinct from the MSSM. 
Nevertheless, let us recall how we can increase the gamma-line signal with features similar to the MSSM. 
In a generic neutralino case, one can look for a point in parameter space where a resonance is hit, 
with m^o — Mz/^ ~ 45 GeV, not forgetting the fact that the ^Z^ channel would require special care 
since we would produce the Z^ boson at threshold. Second a mostly wino, higgsino or a mixture of a 
bino with one of these states would result in a signal with better prospects, like in the MSSM. However 
such type of neutralinos must be quite heavy, especially for wino- or Higgsino-like neutralinos. Indeed, 
if one wants to reproduce the correct CDM abundance, typically the masses must be of order the TeV 
scal^. Motivated by results indicative of light dark matter, we focused on the case of a NMSSM light 
neutralino, which can have very different properties than the one of the MSSM. It is very difficult in the 
MSSM to have a very light neutralino which do not overdose the Universe and at the same time fulfills 
the collider constraints [71|8l[T3l[T3]. Reconciliating these two criteria can be provided by the NMSSM 
more easily. Indeed the lightest of the two CP-odd Higgs bosons can be much lighter than the single 
CP-odd Higgs boson of the MSSM without violating collider constraints. Another feature that cannot 
be reproduced in the MSSM is obviously the one where the LSP is mostly singlino. We have seen that 
on general grounds the signal is very low for such a Xi- Nevertheless, just like with a boson we 
can enhance the signal by hitting a resonance, this time with one pseudoscalar Higgs. This mechanism 
has already been pointed out in |16| . and is also invoked to give the right relic density when the LSP 
possesses an important singlino component [71|30]. We therefore investigated to which extent the signal 
can be enhanced for such a mechanism with neutralinos lighter than 15 GeV, and if the present limits 
can be used to further constrain the parameter space found in [7J. Such light neutralinos would evade 
the present constraints from the published FERMI data on gamma-lines since their threshold is around 30 
GeV pT|. However Vertongen and Weniger [39j have derived a broader energy range which starts from 
1 GeV and extends up to 300 GeV, see Table. 3 of As a consequence we will only discuss how to 

''Scenarios with heavy MSSM wino-hke or Higgsino-Uke LSP also require special care for the gamma-line computation, 
in particular the treatment of the Sommerfeld singularity, see the pioneer work of ^40J 
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increase the 77 signal since the ^Z^ channel is closed in this range, and we will compare our predicted 
rates with the limits given in |39) . 

5.1 Light singlet pseudoscalar Higgs resonance 

Light neutralinos in the NMSSM have been studied for example in [71 l41lH2] . A singlino-like x? emerges 
when K ^ 1 and to ensure at the same time a light singlet pseudoscalar the parameter should 
be small either (in the decoupling limit ~ —"iKA^s). Light singlino-like LSP are found mainly 

in two regions of parameter space : one at small A with very small and another at large A with 
slightly larger allowed, see }41| . To hit a light singlet pseudoscalar Higgs resonance we obviously need 
m^o ~ This can be achieved with a pure singlino-like neutralino but, in practice, neutralinos 

fulfilling this criteria and accounting for the present cosmic abundance of dark matter measured by WMAP 
are mixed states, like a bino-singlino or bino-higgsino like neutralino. In any case the neutralino obviously 
needs some substantial component related to the Higgs sector to couple to the pseudoscalars. Such light 
pseudoscalars and neutralinos can be motivated in models were an approximate R-symmetry or 

Peccei-Quinn symmetry holds (see [41j and references therein) and therefore their "lightness" is 
more "natural" than in the MSSM. To see in which extent the signal can be boosted by this mechanism 
we performed a scan over the parameter spaced spanned by the 11 parameters 

Ml, M2, ^, A, K, Ax,A^, At,mj, nig 

as in [7] and computed the rate (av) ^.^.^-y-y as only this channel is open. The parameters m^- and niq 
are common masses to the sleptons and squarks respectively and are taken in the same range as in [7J. 
However to find these points we did not apply any constraints, we just tried to find points in parameter 
space which where minimizing the difference 

2m^a — niA, 

AM = — ^ (5.1) 

ruA^ 

such that we observed an enhanced signal, and with m^o < 15.5 GeV. No width to the pseudoscalar Ai 
propagator was introduced during the calculation, since we did not encounter any numerical problems. 
This is due to the fact that the width of the pseudoscalar is extremely narrow (numerically 10^'' < < 
10~^). Anyway, as a check, we added a width to the propagator of the lightest pseudoscalar for the points 
giving the highest rates, i.e, very close to the resonance, and we observed a maximum decrease of ~ 20% 
of the signal, so these points would still be excluded. As expected, the effect of the width was less and 
less important (reaching the percent level) the more we were going away from the resonance, where more 
points are allowed. The effects of the width for these allowed points is therefore lower than the present 
experimental accuracy. The result of this scan is displayed in Fig. [2l Red squared points are scenarios 
where the LSP is more than 90% bino, blue diamond-shaped points more than 90% singlino and the 
green stars are a mixture of both. We displayed only the region where |AM| < 1%. We can see that for 
this mechanism to boost significantly the signal a high degree of fine-tuning of AM is required. Very few 
bino-like neutralinos scenarios were found and give a signal lower than 10~^^ cm^ s~^, showing that this 
mechanism is specific to the NMSSM. Most of the points possess also a significant Higgsino component 
which ranges from the percent level to approximately 15%. Notice also that as soon as |AM| > 0.2% the 
rate significantly decreases and more points were produced with rriAi — ^"^x?' observe that in the 
most extreme case we can reach a signal as large as 10~^^ cm^ s~^, for a mixed Xii which is 88% singlino 
and 9% Higgsino while AM ~ 0.001%. However recall that we did not apply any constraining criteria for 
this scan. In Fig. [3] we display the result of this scan with respect to the limits given in |39) . This time 
we picked up points with |AM| < 15%, this selected the mass window 4 < m^o < 15.5 GeV. We see that 



the highly fine-tuned points (with |AM| < 0.2% and (av)^^ > 10~ cm'^/s) would be excluded by the 
spectral lines searches. However we found that most of the points are just below the FERMI sensitivity, 
especially concerning the galactic center observation. Therefore an increase by at least one order of 
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AM[%] 



Figure 2: Gamma-ray line rate with respect to the mass difference AM = (2m^o — rnAi)/'mAi for several 
composition of the neutralino (see text for details). 




mJGeV] 



Figure 3: Gamma-ray line rate with respect to m^. The limits of f&yij are drawn for the halo (dashed 
dotted line) and the galactic center (solid line). The conventions for the neutralino composition follow 
those of Fig. \^ 
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magnitude would much more constrain the s-channel pseudoscalar resonance, if no signal is reported. It 
is also worth to note that the highly fine-tuned scenarios tend also to give a very low relic densitjj^l, hence 
giving up the possibility of the LSP to account for the actual amount of dark matter. Moreover, as we 
will discuss in the next section, the gamma-rays from dwarf spheroidal galaxies (dSph) [431144] give more 
powerful constraints. 

5.2 Gamma-ray lines as a further constraint on the NMSSM parameter space. 

Interestingly the studies led by [71ll5j found that light NMSSM neutralinos that are compatible with 
CoGeNT and satisfy all the constraints they applied (colliders, relic density of DM, flavour observables) 
are accompanied by a light pseudoscalar and/or scalar singlet. Generally the aforementioned studies 
showed that the highest elastic scattering cross-sections are reached for low values of and can satisfy 
the WMAP constraint when mui — 2?7i^o, because in this case the pair annihilation of neutralinos is 
enhanced by a s-channel resonance. Obviously the gamma-ray lines cannot constrain such scenarios as 
this channel is suppressed for low-velocities. Scenarios that could be constrained are therefore rather 
those where m^o ~ 

The authors of [43 j computed the secondary gamma rays produced in dSph from the pair annihilation 
of DM particles into quarks and/or taus which subsequently hadronise and decays into pions, finally 
decaying into photons. The authors then compared their findings with the FERMI 95 % limits on gamma- 
ray emission from dSph |46) . The authors provided us with 14 points of their MCMC scan giving a large 
pair annihilation cross section but safe with respect to dSph limits and direct detection searches. The 
points are sampled in each bins of m^o between 1 and 15 GeV. We then used these input parameters 
to evaluate the rate of the gamma-lines to see if the corresponding limits derived in |39| from the FERMI 
data on spectral lines [17J could further constrain these scenarios. The mass spectrum of these parameter 
points fall into the range where m^i^ is not far from 2m^o . However models for which the mass difference 
between the LSP and the pseudoscalar is highly fine-tuned are already excluded by the Draco limits [43j . 



m-r.0 
Al 

[GeV] 


[cm^ s^^] 


(av).,., X lO^'J 
[cm'^ s""*"] 


0.976 


0.209 


0.00008 


2.409 


0.297 


0.00267 


3.342 


0.345 


0.00345 


4.885 


3.298 


0.00262 


5.626 


5.389 


0.00410 


6.551 


3.547 


0.00427 


7.101 


2.425 


0.00664 


8.513 


2.161 


0.00220 


9.274 


2.497 


0.00655 


10.27 


2.323 


0.01881 


11.50 


2.575 


0.02456 


12.74 


3.224 


0.02003 


13.51 


9.571 


0.17487 


14.48 


148.4 


2.87500 



Table 5: Comparison between dark matter annihilation into quarks and/or taus (values are taken from f^) 
and the loop-induced one into photons in the NMSSM for each of the bins between 1 and 14 GeV. 



^Nevertheless we found some points satisfying this criteria. 
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We therefore do not expect very bright signals. We coUected in Table[5]the total annihilation cross section 
of neutralinos into quarks and/or taus, taken from [43], with our results concerning the annihilation into 
two photons, since here again the range of m^o studied is too low for the production of the ^Z^ final 
state. As can be seen from Tab. O the gamma-line rate is at most four orders of magnitude lower 
than the rate coming from other indirect search channels. Of course since the direct annihilation into 
photons is loop-suppressed we cannot expect that it can significantly contribute to the total annihilation 
cross section producing gamma rays. For the computation of gamma ray fluxes emerging from dark 
matter annihilation in dwarf spheroidal galaxies, the "direct" annihilation into photons can then safely 
be neglected. If we now compare the results in Tab. [5] with the limits of [39], which are at best of order 
10~^^cm^/s for a very light neutralino (less than approximately 5 GeV), none of these points is further 
constrained. Most of the models exhibits a spectrum with a mass difference AM too high to enhance 
significantly the signal, and relatively low mass sleptons. The point with m^o = 14.48 GeV gives the 
best signal because AM ~ 0.5% with a LSP at 92% singlino and 8% Higgsino. A sensitivity increased 
by at least two order of magnitude would be needed to exclude this point with respect to the constraints 
from gamma-lines searches. In consequence we conclude that the present limits on gamma ray lines can 
constrain models with a very low mass difference AM but they are not competing with the gamma ray 
searches from dwarf spheroidal galaxies anyway. This can be easily understood by the fact that photons 
produced from processes like XiXi ~^ ViItt^ are leading-order dominated and experiments have not yet 
reached a sensitivity down to loop effects, like the spectral lines. 

Conclusion 

The mono-energetic gamma ray line signal has spectacular features : a clear "smoking-gun" signature 
and a direct relation to the mass of the dark matter particle. Moreover it do not suffer from astrophys- 
ical uncertainties and depend only on the assumption of the dark matter halo. However discriminating 
it from the overwhelming astrophysical background (supernovae, pulsars, cosmic-rays...) is experimen- 
tally extremely challenging and requires a fine energy resolution. Furthermore the inclusive channels 
XX ~^ qq{WW^ ZZ^ produce a featureless continuous spectrum of gammas coming from the decay or 
hadronisation of the final state, only cut-off at a maximum energy corresponding to the mass of the DM 
particle. This contribution is stronger than the monochromatic photons, since it is a tree-level dominated 
process. In this paper we evaluated numerically in the NMSSM framework the process XiX? ~^ 11 a-iid 
provided for the first time result for XiXi ~^ iZ^ . This is the first implementation of the NMSSM frame- 
work in SloopS and paves the way for a future renormalisation of each sector, like what was done for 
the MSSM |25|I26|. We performed the calculation at f = and checked the result with respect to gauge 
invariance thanks to an NMSSM extension of the non-linear gauge- fixing used in [25]. As a further check 
we compared the NMSSM result in the "MSSM-limit" and the pure MSSM one and found an excellent 
agreement. As a first application we computed the ^^{Z^) signal for several benchmarks found in the 
literature [23p24) . We compared with an "equivalent" MSSM spectrum and the limits given by FERMI |17| 
and in [39]. We found that unless the lightest neutralino has a significant singlino component, the two 
models produce very similar rates and discriminating between them requires a high level of experimental 
precision. This is particular true if in both cases the LSP is bino-like. However, if the neutralino is to 
be the dark matter candidate and light, discriminating the NMSSM from the MSSM should be easier 
since light neutralinos and Higgs scalars are more "natural" in the NMSSM, thanks to an approximate 
Peccei-Quinn symmetry. A low mass WIMP is therefore easier to accommodate in the NMSSM than the 
MSSM. In the NMSSM the Higgs sector is more peculiar and its related collider observables are expected 
to give also more spectacular signatures, hence are more prone to unveil the singlet component of the 
Higgs potential [15). Therefore dark matter related observables in the NMSSM are also distinct when 
the singlet component is at play both in the neutralino and Higgs sector. Furthermore, as far as the 
annihilation XiXi ~^ ll{^^) is concerned, this observable can only be sensitive to the pseudoscalar Higgs 
states, since at f = scalar couplings are suppressed. In this particular case, such regions of parameter 
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space are only reachable in the NMSSM. A noticeable example of such a mechanism is self-annihilation 
of dark matter particles through a s-channel resonance with a light pseudoscalar. We investigated such 
a mechanism with a light x? for the 77 final state and observed that the signal can be very bright when 
extremely close to the resonance. A singlino-like LSP not fulfilling this feature would on the contrary give 
a very low rate. These kind of highly fined-tuned scenarios would in fact be excluded taking into account 
the limits derived by the authors of [39j- However, if this mechanism is actually at play for the indirect 
detection of primaries photons, it is also the case for both the diffuse photon spectrum and relic density 
calculations. Indeed, if ones wants to match the CDM abundance, this fixes the pair annihilation rate 
to be of order the canonical value of {av) ~ 3 x 10^^'^ cm^/s. Interestingly, this value is in the ballpark 
of the FERMI sensitivity. Unfortunately, as the gamma lines is a loop suppressed process, the interesting 
sensitivity is one or two order of magnitude less, in the most favourable scenarios, like the pseudoscalar 
resonance. Outside these regions, limits on the monochromatic gamma lines are not very constraining, in 
particular, on the astrophysical level, with respect to the ones from dwarf spheroidal galaxies. Further- 
more, even concerning the pseudoscalar resonance, the dSph put also some more stringent constraints 
on this mechanism, since it is also at play there. The sensitivity of experiments concerning the spectral 
lines still need to be improved by several orders of magnitude to be compelling. The FERMI collaboration 
claims that their limits on the spectral lines will be significantly improved by the time the mission is 
finished and we can expect a possible discovery or, in case of a null result, that more featureless regions 
of the NMSSM could be excluded. 
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Appendix : NMSSM BRST transformations for the scalar fields 

To generate the ghost Lagrangian through the BRST transformations starting from the non-linear gauge- 
fixing in Eq. ()2.ip we need, in addition to the BRST transformations for the gauge fields which can be 
found in |35| . the BRST transformations of the scalar fields from the gauge transformation of the Higgs 
doublets. We can parametrise the fields of the Higgs sector in the following way 



The following unitary rotations matrices turn the gauge eigenstates to the physical ones 

CP - 

h 




(A.2) 
(A.3) 



'Pi 

(^2 I = I Z'^,C(, Z^,c^ -s^ ] \ A2 ] , (A.4) 





fHi\ 


)l 











^21'5/3 


C/3 ^ 




( Ai\ 


Zficp 






II 




ya 


ya 
^22 









(A.5) 
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where we have first rotated with an angle /3 in the pseudoscalar sector to single out the goldstone mode 
from the physical CP-odd states in Eq. ()A.4p [15j. As the BRST transformations are closely linked 
to gauge transformations we have for the singlet components, 



Sbrs^s = (A. 6) 



Provided the rotations matrices we have then a relation between the BRST transformations of the physical 
eigenstates once the ones for the gauge eigenstates are known, 



With 




-Sj3 Cp ) \ 5BYis4>2 

yh yh yh 

^11 ^21 ^31 

7h lyh lyh 



\ / ^BRS<; 


M 






/ V 






(A.7) 



^12 ^22 ^32 I I '5brs02 | (A.8) 
yh yh yh 
^13 ^23 ^33 

^f2 I I '^BRSV'2 I (A.9) 
C/3 



^ie [c^ - fi^") [c^G± - SfiH^\ (A.IO) 

5BKS(t>t = [^2 + {Zl^lHi + Z2^2^2 + ^2^3^3) ± i(C/3(Z?i^l + Zfi^a) " SpG^)] 



^ie ( - £^^^ ) [s^G± + c;3^±] (A.ll) 



52 _^2 

'2SyjCi 

5brs</'? = +|[c;3(G-c+-G+c-)-s;3(^-c+-F+c-)] 

+-^c^[s^(Zfi^i + Zfi^a) + c;3G°] (A.12) 

8Bns<t>l = +'-^[sfi{G-c+ -G+c-) + cp{H-c+ -H+c-)] 

-c^[cp{Zt^Ar + Z^^A2) - spG^] (A.13) 



6bksV>i = +^[cfiiG-c+ + G+c-)-sp{H-c+ + H+C-)] 



-c^[vi + Z'l^Hi + Zf2i^2 + ^{3^3] (A.14) 



w 



'^brsV'2 = -\[sp{G-c+ + G+c-) + cp{H-c+ + H+C-)] 

+ -^C^[V2 + Z2^//l + Z2^2^2 + ^2^3^3] (A.15) 
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